NASA CR-152020 





(NASA-Cri- 152C20) COMtUIEB TBANSFORMATICN Of N77-28d66 

PAHIIAL DIffEEEMiAL tgUATICNS INTO ANY 
COOHEINAXE SYSTEM (Aercnduticai Pesearch 

Associates oi Ptircetcii) 41 p HC AC3/MF A01 Unclas 

CSCL 12A G3/64 37223 






AU'’- 1S77 


RECE'VED 

^ t NASA sn FAWUIY 
V INPITT BRANCH 


H 


A. R. A. P 


V.. 


s 






AERONAUTICAL RESEARCH ASSOCIATES of PRINCETON, INC 



NASA CR-152020 


A.R.A.P. REPORT NO. 309 


COMPUTER TRANSFORMATION 
OF PARTIAL DIFFERENTIAL EQUATIONS 
INTO ANY COORDINATE SYSTEM 


ROGER D. SULLIVAN 


JUNE 1977 


PREPARED UNDER CONTRACT NO. NAS2-8014 


BY 

AERONAUTICAL RESEARCH ASSOCIATES OF PRINCETON. INC. 
50 WASHINGTON ROAD. PRINCETON. NEW JERSEY 08540 

FOR 

AMES RESEARCH CENTER 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 



ABSTRACT 


Tensors provide a compact way of writing partial differ 
ential equations in a form valid in all coordinate systems 
However, in order to find solutions o'f the equations with ’ 
their boundary conditions they must be expressed in terms of 
the coordinate system under consideration. For complicated 
coordinates the expressions are very long and are tedious to 
calculate. The process of arriving at them from the tensor 
formulation has been automated by a software system, TENSR, 
described in this paper. An allied system that analyzes the 
resulting expressions term by term and drops those that are 
negligible is also described. 
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1 . INTRODUCTION 


It is natural when computing the flowfield about a cone 
to use a conical coordinate system. The description of the 
body surface in that system is simply the statement that one 
of the coordinates is equal to a constant, for example 

6 = 0c 

This makes the application of the boundary conditions to the 
equations of motion straightforward at the cost of compli- 
cating the equations themselves. The extra complication is 
worthwhile compared to the computational problems that arise 
when boundaries cut across coordinate surfaces. 

Such a scheme can be regarded in two ways. In one view 
the conical coordinate system is an entity in Euclidean space 
which is more convenient than the Cartesian system for the 
problem at hand. In another vieM-it is a mapping of the body 
and the flowfield around it into a rectangular grid. To 
illustrate with a two-dimensional example, a polar coordinate 
system can be viewed as a way of locating points in physical 
space, as in this sketch: 
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or as a mapping from that space to this: 



The mathematics is the same in both cases, but the mapping 
interpretation makes it manifest that one gets the advantages 
of finite-difference methods for a rectangular grid along 
with simple boundary conditions. 

Coordinate systems or mappings can also be devised so 
that points are unevenly spaced in physical space (so as to 
be concentrated near a body, for example) but mapped to evenly 
spaced points in the rectangular grid, again gaining an advan- 
tage for finite-difference computations. 

The gains, however, are generally at the cost of more 
complicated equations to be solved. Nevertheless, the advan- 
tages of such a procedure for computation are so great that 
Thompson, Shanks, and Walker (Reference 1) use a time-dependent 
version that keeps a free surface, in the problem they consider, 
mapped into a fixed coordinate line. Thus at each time step 
of their calculation they recompute the mapping as an aid in 
computing the flow. . 

So far, only transformation of the independent variables 
of a problem has been mentioned. It is often desirable to 
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transform the dependent variables, in the sense that the repre- 
sentation of vectors changes from components aligned with the 
original coordinates to components aligned with the new coordi- 
nates. For example, in transforming from Cartesian to polar 
coordinates in the plane, the representation of velocities 
would change from horizontal and vertical components to radial 
and tangential components. This also is ordinarily advantageous 
for the expression of boundary conditions. Its effect on the 
complexity of the equations is varied. There is a gain in 
unifying the representation of dependent and independent vari- 
ables, but on the other hand new terms arise. The prototype 
of such terms is the centripetal acceleration term which turns 
up when radial and tangential velocity components are used. > 

At first glance it would seem that the transformation of 
the dependent variables as well as the independent variables 
would greatly complicate the task. However, the techniques of 
tensor analysis accomplish both transformations in one step 
with great economy. On the other hand, the compact tensor 
notation can disguise very long expressions which must be ex- 
panded before much can be accomplished numerically. (For tensor 
expressions which Involve known quantities there is no problem 
in writing programs to calculate them. It is expressions 
involving unknowns that are principally considered here.) 

As the body shapes considered grow more complicated, the 
coordinate systems adapted to them grow more complicated, the 
equations of motion expressed in terms of them grow longer, and 
the likelihood of getting them written down correctly decreases 
if the operation is done by hand. But the process at this 
stage is purely mechanical and so can be done by computer. The 
software system TENSE accomplishes Just that. 


The original version of TENSR, though, was produced to 
answer a different need. Since 1967 A.R.A.P. has been engaged 
in second-order closure techniques for calculating turbulent 
flows, and since 1968 the modeling has been done in the frame- 
work of tensor analysis in order to eliminate the bias that 
can easily be Introduced when working in a fixed coordinate 
system. The principal unknown in our system of equations for 
turbulent flow, besides the variables describing the mean flow, 
is a tensor of second rank, namely the correlation of component 
of the velocity fluctuations, • '^he governing equa- 

tions are quite complex even in Cartesian coordinates. We 
achieved good results solving these equations for incompress- 
ible flows (Reference 2) but not without considerable problems 
in accurately developing computer programs for the purpose. 

Therefore, when ^^^e approached the task of solving the 
equations of turbulent compressible flows — for which fluctu- 
ations of density, temperature, viscosity, and so forth, as 
well as of velocity, must be taken into account — we decided 
to automate the process even though only Cartesian coordinates 
were to be used, since the study involved the boundary layer 
on a flat plate (Reference 3). Thus, TENSR was born able to 
handle only Cartesian cooi’dlnates . 

The output of the original TENSR was edited by hand to 
remove terms that are small in the context of boundary-layer 
theory. After TENSR was extended to be able to handle arbi- 
trary systems, more complicated coordinates were processed 
and the output grew so much that it was desirable to make the 
dropping of terms automatic. For this purpose, a follow-on 
system, TATTR (for Term ATTRlbutes), was written. The use of 
this program not only eased the burden of dropping terms but 
also forced the process to be more uniform and logical than 
had sometimes been the case in our work. 


Another software system, DIFPR, uses the edited output 
of TENSR or of TATTR to create Fortran statements which become 
a main part of programs that solve the equations by finite- 
difference techniques. 

This paper reports the operation and characteristics of 
TENSR and TATTR. Since the present version of DIFPR is tied 
to a specific finite-difference scheme, it is not further 
described here. 

In the next section the notation used to handle tensors 
on the computer is introduced in the course of a brief lesson 
on tensor concepts. Sections 3 and 4 show some examples of 
what TENSR and TATTR do. These are followed by a section on 
possible future developments and on related efforts by others. 
Some characteristics of TENSR that may interest tensor experts 
are considered in Section 6. Finally there is a brief bibli- 
ography. 


2. TENSORS: A QUICK LESSON 


One view of tensors, sufficient for the understanding of 
TENSR, is given here. Readers already familiar with the sub- 
ject are advised to skim the section to gain familiarity with 
the notation we have devised to communicate tensor expressions 
with computers. As is explaiend in Section 6, TENSR is not 
restricted in the type of space it deals with (except for the 
upper limit of ten dimensions) but applications such as those 
mentioned in Section 1 are met with in an ordinary three- 
dimensional Euclidean space so this lesson is confined to 
that case or a two-dimensional subset of it. The treatment 
here is expository; proofs may be found in the works cited in 
the bibliography. 


In the powerful index notation used in tensor analysis 
the equation 


Cl = 


( 1 ) 


stands (in three dimensions) for three equations each with 
three terms on the right : 

Cl = A^^B^ + A^^B^ + A^^B^ 


°2 = 


( 2 ) 


Cg = A^gB^ + A^^Ba ^^3^3 


The rules governing symbolic indices (that, is , non-numeric 
Indices) are these : 

• A symbolic index appearing alone in a term, known as 
a free index, takes on the range of values in separate 
Instances . 


• A symbolic index appearing twice in a term, known as 
a dummy index, implies a summation over the range of 
values (the Einstein summation convention). 

• A symbolic index doesn’t appear three or more times 
in a single term without a special Indication of the 
meaning to be assigned to the expression. 

Thus, the usual apparatus that would accompany the con- 
densation of equations (2), for example: 

3 

Oj^ = y (1 = 1.2,3) 

n=l 

is eliminated, leaving just equation (1). 

Some additional rules, important in the general case, but 
unnecessary and therefore often disregarded when only Cartesian 
coordinates are considered, are: 

• Dummy indices in a term appear once up (as a super- 
script) and once down (as a subscript or in a 
derivative) for each pair, 

• Free indices appear in the same position in each 
term in an equation. 

The comparison of equations (l) and (2) gives only a hint 
of the economy of this notation* In three dimensions an equa- 
tion with three free indices, e.g., equation (7) below, stands 
for 27 equations (not necessarily all independent) and a term 
with four pairs of dummy indices, e.g., in equation (8) below, 
stands for the sum of 81 terms. 

In devising a notation to use on the computer, we decided 
to adopt conventions like those of Fortran in many respects. 
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but the Fortran subscript representation was discarded since 
it lacked a means of representing superscripts and since it is 
unwieldly when expressions teem with indices. Instead, we 
decided that indices would consist of only one letter or digit 
and that superscripts would be preceded by " and subscripts 
by ' . Thus, equation (1) appears as 

C‘I = A"N‘I*B»N 

Since variables are not restricted to one character, this m.ight 
appear as 

CRS'I = AAA1’'N'I*BS»N 

Several superscripts or subscripts need not have separate " 
or ’ . Thus, A^J is represented by A"TJ . 

Let a Cartesian coordinate system be denoted by y^ and 
some other system by . Functions relating one to the other 

are assumed to exist and in most applications are known in 
advance. For an example in tvfo dimensions, polar coordinates 
are related to Cartesian coordinates by 

11 2 
y-L = x-^ cos x“^ 

y2 = xl sin x^ 

Actually, the x^ are seldom used in this context; it is more 
convenient to use r and (j> : 

y^ = r cos cj) 

y2 = p (j) . 

But in speaking of the coordinates collectively, or especially 
in speaking of coordinate systems generally, the notation x^ 
is used. 



Scalar fields, that Is, scalars which are functions of 
position, such as density, are tensors of rank zero. The 
gradient of such a field is a tensor of rank one, that is, a 
vector. Vectors, like scalars, are invariant with respect to 
coordinate transformations. However, the representation of 
a vector in components does vary with the coordinate system. 

A representation in components which transform as do the 
partial derivatives of a scalar field is called a covariant 
representation of the vector and the indices appear as sub- 
scripts. For example, we might set 

C, = ^ ^ . 

(The colon is used to denote partial differentiation both in 
normal notation and in TENSR notation: C’l = PSI:I .) a 

representation in components which transform as do differential 
displacement components is called a contravarlant representa- 
tion of the vector and the Indices appear as superscripts. 

For example 




dt 

Covariant representations can be calculated from contra- 
varlant representations and vice versa by the formulas 


Ai 

= 

or 

A»I 

aJ 


or 





This is called lowering or raising the index. In place of the 
phrase covariant (contravarlant) representation of a vector, 
one often speaks concisely of a covariant (contravarlant) 
vector. Also, an index is often called Covariant or contra- 
varlant according to its position as a subscript or superscript. 

10 - ; , . 


As a mnemonic, note that sub- and op- are short, whereas 
super- and contra- are long. It should be clear that covariant 
representations of vectors are not necessarily partial deriv- 
atives of scalars nor are contravarlant representations of 
vectors necessarily differentials or velocities. 

Higher rank tensors may have representations which are 
covariant (all subscripts), contravarlant (all superscripts), 
and mixed (some of each) and the indices may be raised and 


lowered individually. 

For example : 



= D;^^rn 

® ^nj 

or 

T"I ’ J 

= @’'IN^fT'NJ 

m 3 


or 

Tvr'j 

= @"JL*T'IL 

ijii J 


or 

H 

II 

@''JL#@"IN«T»NL 


Equation (3) can be looked at as raising the first index of 
glj or lowering the second of g^J ; in either case, we 
conclude that 

or @"I’J = &"i»j 

That is, the value of the mixed representation of the metric 
tensor is the Kronecker delta in any coordinate system. This 
contrasts with the fact that the value of g. . is 6 . and 
the value of g^^ is 6^^ In Cartesian coordinates only. 

If a superscript and a subscript of a tensor are made the 
same, thus implying summation, a new tensor is formed of rank 
two less than the original — a process called contraction. 
Thus, the two vectors 


_ - mn 

ni 

or 

A» I 

II 

»NI 

rpn 

^ in 

or 

. B’l 

= T”N 

’IN 
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are different (in general) contractions of or T"I*JK 

Some tensors are the product of tensors. A tensor of rank 
multiplied by one of rank q yields a tensor of rank p + q 
The contraction of the product of two vectors, or 

A"N»B'N Is the scalar product of the vectors. 


P 


In addition to covariant and contravarlant components 
of tensors there are also physical components, that Is, com- 
ponents In the elementary sense. Velocity In polar coordinates 
provides an Instructive example. The physical components are 
denoted by u for radial and v for tangential, the covariant 
and contravarlant components by Uj_ and u^ . in this system 
the radial component Is the same In all three versions, 
covariant, contravarlant, or physical: 


= U^ = U 

On the other hand, the three types of tangential component are 
quite different. The contravarlant component can be described 
as an angular velocity, since 



whereas the covariant component Is more a circulation, since 

''2 = ( 6 ) 

See References 4 and 5, In particular, for more on physical 
components. 


The representation in polar coordinates of the velocity 
of a body moving In a circle at constant speed is unchanging 
That is, the radial component Is constant (at zero) and the 
tangential component Is constant. Nevertheless, as is well 


known, the body is accelerating toward the center and the equa- 
tions of motion in polar coordinates represent the situation 
by the centripetal acceleration term. The generalization of 
this term to situations Involving arbitrary changes in tensors 
using arbitrary coordinate systems is embodied in Chrlstoffel 
symbols. The Chrlstoffel symbol of the first kind is given by 

^ijk ~ ^ (Slj;k ®ik:j “ Sjk:i) 
or ( 7 ) 

$’IJK = (@’IJ:K + @*IK:J - @'JK:l)/2 

(The first choice for TENSR notation for the Chrlstoffel symbols 
was (f: , but A.R.A.P. * s line printer didn’t have it available 
so $ was chosen as the next best thing. Since the literature 
is not in agreement on notation for them and since the ordering 
of indices that I find convenient — that is, such that is 

symmetric with respect to j and k — Is unorthodox, I have 
adopted the practice of using $ in ordinary as well as TENSR 
notation.) The Chrlstoffel symbol of the second kind is 

^^jk S^^$njk or = @"IN»$»NJK 

For our purposes, the first kind has no function except to allow 
the calculation of the second kind. Neither kifid of Chrlstoffel 
symbol is a tensor. 

Covariant differentiation, denoted by a comma preceding an 
index, is an operation that produces tensors from tensors. When 
applied to a contravariant vector (to slip into that terminology) 
its rule for calculation is 
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Neither term on the right is a tensor, but the sum is. The 
covariant derivative of a covariant vector is 


^i,j = 


The generalization to tensors of higher rank is straightforward. 
For example : 


^ k,r ^ k:r 


*i A*^^ + 

^ ar^ k ■ 




ia *a .ij 
k “ ^ kr^ a 


The rule is not difficult. Each index in the original tensor 
is replaced by a dummy in one of the extra terms. The rest of 
the index juggling is automatic by the rules for indices. The 
only aspect that is not automatic is the sign thS-t accompanies 
each extra term. (The user of TENSR need not bother with any 
of this — TENSR does it all automatically.) 


Since a scalar has no Indices, its co variant derivative 
is simply its partial derivative: 

''.I = ''ti 

This agrees with the previous assertion that the gradient of 
a scalar is a tensor. In Cartesian coordinates the Christoffel 
symbols are all zero; therefore, partial derivatives in 
Cartesian coordinates are also covariant derivatives. 


It is convenient to have a notation for what might be 
called a contravariant derivative, that is, a covariant 
derivative with the index raised. Again, the literature is 
not consistent; TENSR uses an exclamation point: 


PSI!K = (§"KM*PSI,M 



This brief lesson misses much of the splendor and impor- 
tance of tensors. The reader is urged to consult some of the 
works cited in the bibliography for a better view. It should 
be noted that many of those books start by considering abstract 
spaces for which a metric does not exist. In that case there 
can be no raising or lowering of Indices, and so covariant 
and contravarlant vectors are truly different types of entitles, 
unrelated to each other. Once spaces with metrics are con- 
sidered, however, the viewpoint becomes closer to that presented 
here. 



3. SOME TENSR SAMPLES 


In ordinary applications, TENSR accepts information con- 
cerning the metric tensor and Chrlstoffel symbols for a coordi- 
nate system, along with symbol substitutions to be made, and 
one or more equations in general tensor form. It proceeds to 
expand (internally) covariant derivatives into the sum of 
partial derivatives and Christoffel-symbol terms, to expand 
Implied summations, to drop terms which are zero, to make sub- 
stitutions specified, and to print the results and store them 
on disk. 

For a short example, consider the transformation to polar 
coordinates of the Navier-Stokes equations for the steady flow 
of a fluid of constant density. Figure 1 shows the printout 
from such a run. The first seven lines are exact images, of 
input cards. The symbol # is used to indicate comments; any- 
thing to the right of it is Ignored by TENSR. The first card 
sets the range of dummy index summations as 1 and 2, and hence 
establishes that a two-dimensional space is to be considered. 

The next two cards define the metric; since the metric tensor 
is symmetric only three of the four elements need be supplied 
in each case. A blank followed by the separator ; suffices 
to indicate a zero element. Fortran conventions concerning the 
order of operations are observed; thus 1/R/R = (1/R)/R stands 
for r“2 , The fourth card supplies values for the Chrlstoffel 
symbol of the second kind. Again because of symmetry, only six 
of the eight values are specified. Symbol substitutions are 
designated by the next card* In this case expressions involving 
physical components are to be substituted for covariant and 
contravariant components [compare equations (4)-(6)]. The 
second DEFINE INDEX card identifies K as a free index; in the 
absence of other indication the range of the free index is the 
same as the range of dummy indices. Finally, the equation to 


DEFINE INDEX 1 , 2 


DEFINE S" l; ; 1/R/R 

DEFINE a' l; ; r*r 

DEFINE ; ; -r; ; 1/r; ; 


^ POLAR 


POLAR 


# POLAR 


DEFINE SYMBOL U"1 = U; U”2 = V/R; U'1 = 


u; U'2 = RAV # POLAR 


DEFINE INDEX K 


U"L*U'k,L = -P,< t NU*a"LN*u’K,LN # NA V I E R -S T OKE S 

U"L*U’K/L = -P/K + NU*a"LN*U *K ^LN 

U*(U):1 + V/R*(U):2 - v/R*1/R*R*v = -(P ).1 + NU*( (U) :1 ) :t 

+ NIJ*1/r/r*((u):2):2 - NU* 1 /R / R * ( 1/ r* r* v) ; 2 - NU*1 /R/R*1 /R* (r*v) : 2 

- NU*1 /R/ R * 1 / R*R*U + NU*1 / R /R *R* ( LI) : 1 

U*(R*V):1 + V/R*(R*V);2 - IJ*1/R*r*v + V / R* R* U = ( P ) : 2 

+ Nl)*C(R*V):1):1 NU*1/R/R*((R*y):2):2 - NU* ( 1 /R* R*V ) : i 

+ NU* 1 / R/ R* ( R*U) : 2 - NU*! /R*(R*\/) : 1 + N U* 1 / R / R ★ R ★ ( U ) ; 2 
+ NLJ*1/R*1/R*r*\/ - NU*1/R7R*R*1/R*R*v + NU*i/R/R*R*(R*V):r 

- NU*1 /R/ R*R*1 /R*R*V 


Figure 1 , The Navler—Stokes equations in polar coordinates. 



be expanded appears. Here P stands for the pressure divided 
by the density and NU for the kinematic viscosity. 

After TENSE checks for errors such as unbalanced paren- 
theses, it confirms its verification by reprinting the equation. 
Then the results appear: the radial momentum equation and the 

tangential momentum equation, that is to say the Navier-Stokes 
equations in polar coordinates. There in a nutshell is the 
action of TENSE. 

That TENSE is weak algebraically is evident. Some of 
these terms can be combined and many can be simplified. For 
example, the third term in the radial equation is hardly recog- 
nizable as the centripetal acceleration term since it appears 
as a transcription of the original form, -U'’2*$"2 • 12*U^ 2 , 
without simplification to -V»V/E . In the original appli- 
cations of TENSE, involving Cartesian coordinates, this short- 
coming was very minor and so was allowed to stand. When it 
later became a problem, it was more convenient to deal with it 
in TATTE, the follow-on to TENSE, rather than in TENSE Itself 
( see Section 4 ) , 

To illustrate the functions of the DEFINE cards. Figure 2 
shows the output of a run with all but the first omitted. The 
equation to be expanded is the same but the result is quite 
different. Here, there are none of the replacements and 
dropping of terms that took place in the run illustrated in 
Figure 1. The free index, K , appears as such in each term. 
For each value of K (1 and 2), this represents one of the 
two-dimensional Navier-Stokes equations for arbitrary coordi- 
nates. 

The explosion of terms that can occur begins to show here. 
Applying the rules for covariant derivatives given in Section 2, 
the viscosity term can be written 
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DEFINE INDEX 1 , 2 


U"L*U'K#L = -P/K + NU*a)"LN*U' K/LN H N A V I ER - STOKE S 

U''L’*'U'K/L = -P/K + NU*a"LN*U *K/LN 

U"1*CU'K):1 + U"2*(U'K):2 - U'*1 *$*’1 * K1 *U ' 1 - U*'1 *$*'2 ' K1 *U' 2 

- U"2*$'n ' K2*U ' 1 - U‘’2*$"2'K2*U'2 = -(P):K + NU*3" 1 1 * ( ( U ' K ) ; 1 ) : 1 
+ NUAffl" 1 2* ( (U ' K) : 1 ) : 2 + NU*o)"21 * ( (U 'K ) : 2) : 1 + NU*a"22* ( (U* K) : 2) : 2 

- NU*a" 1 1 * ($"1 'K1 *U ’ 1) : 1 - NLI*a"1 1 * {$"2 ' K1 *U • 2)i: 1 

- NU*a" 12* <$*'1 ' K1 *U ' 1 ) : 2 - NU*a"1 2* ($"2 • K1 *U • 2) : 2 

- NU*B"21 * ($"1 *K2*U’1 ) : 1 - NU*a"21 * ($"2 ' K2*U ' 2) : 1 

- NU*3*’22* ($"1 'K2*U' 1 ) :2 - NU* 9"2 2* ($ ”2 ’ K2*U ' 2 ) : 2 

- NU*a"11*$"1 'K1*(U' 1 ) : 1 - MU*a"11*$"2'Kl*CU'2) :1 

- NU*a*’12*$"1*K2*(UM):1 - NU*a ’*1 2'*$'' 2 ' K2 * ( U * 2 ) : 1 

- NU*a"?1 *$"1 'K1 * <U' 1 ) : 2 - NU*a"21 *$"2 'K1 * (U' 2) :2 

- NU*a’'22*$"1 'K2*(U'1 ) :2 - NU*3 "22*$ " 2 * K2 * ( U ' 2 ) : 2 

+ NU*a" 11 *$"1 'K1 *$"1 • 11 *U • 1 + NU*a"11*$"1 'K1*$"2'11*U'2 
•f NU*a" 1 1*$"2 ' K1 *$"1 ' 21 *U • 1 + NU*a"11*$"2'Kl*$"2'21*U'2 
+ rgU*a"12*S"1’K2*S"1M1*U'1 + NU*a"12*$"1 ' K2*$"2 ' 1 1*U’ 2 
+ NU*a"12*$"2'K2*$’’1'21*UM + [\lU*a"12*$"2'K2*i"2’21*U'2 
+ NU*a"21*$"1'K1*$"1M2*U'1 + NU*a"21*S"1'K1*$"2'12*U'? 

+ NU*a"21*$"2*Kl*$"1 * 22*U ■ 1 + NU *3" 2 1 * $ " 2 * K 1 *S " 2 * 2 2 * U * 2 
+ NU*a"22*$"1 ’K2*$"1 ' 12*U • 1 + NU*a"22*$"1 ' K2*$"2 M 2*U ' 2 
+ NU*a"22*$"2'K2*$"1‘22*U'1 + NU*a"22*$"2'K2*$"2'22*U'? 

Figure 2. The Navier-Stokes equations for general coordinates 
in two dimensions. (Continued on next page) 
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” N U * a " 

- NU*5)"12*$"1 M2* (U»K) : 1 

- NU*a'*21*$*'V21*(u'K) : 1 


- NU*a"11*$"2Ml*(UVK) ;2 

- NU*a''12*$''2 M2*(U'K) :2 
~ NU*3 "21 *$" 2 ' 2 1 * ( L) ' K ) : 2 


- NU*a"22*$"1 *22* (U*K) : 1 - 
+ NU*a"11 "-J *1 1*$" 1 • ki*u m 

+ Nu*a"ii*$»2 . • k2*u' 1 
+ NU*a"12*$"1 M2*S"1 'K1*UM 
+ N'U*$"12*$"2M2*$"1 *K2*UM 
Nu*a"2i*$"i '2i*$"i 'ki*um 

+ NU* a "21 *$"2'21*$"1 ' K2*U M 

* NU*al"22*$"1 • 22*$"1 »K1*U M 
+ NU*a:"22*$"2'22*$"1 'K2*UM 


NU*a"22*$"2 • 22 * CU 'K ) : 2 
+ NU*a"11*$*'Ml1*$"2'K1*U'2 
* NU*a"li*$"2M 1*$"2 'K 2 *U' 2 

+ NU*Qj"12*$"1 M2*$"2‘K1*U'2 
+ NU.*a" 12*$"2 M 2*$"2 * K2* U ' 2 

+ NU*a"21 *$"1 ' 21*$"2 ' K1 *u' 2 

+ NU*a"21 *$"2 '21 *$"2 •K2*U'2 
+ NU*a"22*$"1 •22*$"2'X1*U'2 
+ NU*a"22*$"2 * 22*$ "2 ' K 2*U ' 2 



Figure 2 . 


The Navier-Stokes equations for general coordlnat, 
In two dimensions. (Continued ffom preceding p^g! 


20 




k.*n = ''s'X-ta - 




kn“6:Jl '’®*"^^kn*“en“ct 


!ase ‘tl'ortT *he two-dl.enslonal 

case, two of these terms yield 16 terms each, three yield 8 

o 77 JTl ^ 60 out 

arms n the equation. (The equations of Figure 1 

e\^p7v:r--^r::d":r77^^ °7\7 --ois 

simple polar coordinate sysLm. ) % the77 cird 7thls 

example had been omitted too, the default range for dummies, 

. 2 . and 3 , would have applied. The result would have been 

an equation of 265 terms (252 of them viscosity terms): one 

Of the Navler-Stokes equations In arbitrary coordinates for 

three dimensions. Since this Is one of the simplest of the 

equations used to study fluid flnwc; +-v,/:a -i 

y j,uiQ iiows, the usefulness of TENSR 

can be appreciated. 

Another example Introduces some more features. Figure 3 
s ows a run deriving the axial momentum equation for an axlsym- 
metrlc^jet with no swirl. The cylindrical coordinate system 
given by x-l = R , x"2 = * , x-3 = 2 la used and the 

corresponding physical velocity components are denoted u V 

Plgur 7 between the Input to this run and that of ’ 

Igure 1 are: (1) The range of the dummy Indices Is set by 

e ault. ( 2 ) There are more elements of the metric tensor and 
o the Ohrlstoffel symbol to establish, and there are more 
su stltutlons to be specified. (3) The continuation of the 
preceding line Is Indicated by a character In column 6 as In 
or ran. (l|) The DEFINE NULL card establishes' axlsymmetry by 
Indicating that derivatives with respect to 4. (:2) are zero 
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DEFINE 

a" i; ; ; 

1/r/r; 

; 1 



CYLINDRICAL 

DEFINE 

a » 1 ; ; ; 

r*r; 

; 1 




CYLINDRICAL 

DEFINE 

fi • • • 

^ # 9 9 

* R / 9 

/ 



ft 

CYLINDRICAL 


; 1 / R ; 

9 9 

/ 9 



n 

CYLINDRICAL 


9 9 9 

9 9 

9 



u 

CYLINDRICAL 

DEFINE 

SYMBOL U"1 = 

u ; u 

”2 = 

v/R ; 

u"3 = w ; 

u 

CYLINDRICAL 


UM 

u ; u 

'2 = 

R*v ; 

U»3 = W 

n 

CYLINDRICAL 

DEFINE 

NULL :2 ; U 

"2 ; u 

•2 





DEFINE 

INDEX K = 2/ 

3 







U"L*U'K^L = -P/K + NU*U*K!L/L # N A V I E R-S TOKE S 

U"L*U*K/L = -P/K + NU*U’K!L/L 

0 = 0 —:—-, : 

U*(W) :1 + W*(W) ;3 = -(P) : 3 + NU*( (W) : 1 ) : r + NU* ( (W) :3) ;3 
+ MU*1 /R*(W) ; 1 


Figure 3. The axial momentum equation in cylindrical coordinates. 
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and specifies no swirl by indicating that U”2 and U'2 are 
zero. (This being the case the presence of U"2 and U'2 
on the DEFINE SYMBOL cards is superfluous, but it is convenient 
to have them there for other applications.) (5) Since not all 
three values of K are desired, the DEFINE INDEX card specifies 
are. For the axial ecjuatlon only, K — 3 would appear, 
but the value 2 has been included to show that TENSR handles 
this extreme case, 0=0, correctly. (6) The viscosity term 
has been written with the ’'contravariant derivative" symbol, 

! , instead of explicitly using the metric tensor as in Flgy.r®c 

1 and 2. 


These examples show the general features of TENSR and 
demonstrate that it produces correct results. They only hint 
at its necessity j the modeled eq.uatlons for turbulence in use 
at A.R.A.P. generate thousands of terms for quite simple 
orthogonal coordinate systems (Reference 6). 
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4. SOME TATTR SAMPLES 

The equations of Figure 1 were submitted to TATTR with the 
results shown in Figure 4. In addition to the equations, TATTR 
was given the following Information regarding the dimensions 
(in the sense of dimensional analysis) of the quantities that 
appear as factors in those equations; U and V are velocities, 
P (being pressure divided by density) is a velocity squared; 

NU is a velocity times a length; R is a length; differenti- 
ation with respect to R (:1) is equivalent to dividing by a 
length; and differentiation with respect to $ (:2) is dimen- 
sionless. These were supplied as the exponents to be attached 
to a length and a velocity: 1 and 1 for NU , G and 2 for P , 

and so forth. TATTR has simplified and combined terms, as can 
be seen by comparing Figures 1 and 4, However, the simplifi- 
cation is not complete since TATTR does not recognize that R 
can be taken out of a $ derivative; for example, the seventh 
and eighth terms in the first equation could be combined. 
Nevertheless,' the Improvement is significant. 

TATTR has also printed the exponents giving the dimensions 
of each term next to it. It is seen that both equations are 
dimensionally consistent. This is rather trivial in such a 
simple example, but the ability of TATTR to make this check has 
uncovered errors in long equations that might otherwise have 
gone undetected. 

The dimensions of a term are examples of one type of 
attribute that TATTR can evaluate. Any concept that can be 
quantified by Integers attached to the factors that appear, 
such that the integers add (or subtract) when the factors 
multiply (or divide) to form a term, is valid as an attribute 
for the purposes of TATTR. The type of attribute which motivated 
the design of TATTR is order-of-magnitude information. To 
Illustrate, the axial-momentum equation of Figure 3 was used* 
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LENGTH VELOCITY 


U* (U) : 1 

+ 1/R*V* (U) : 2 
- 1/R*V*V 
= -( P) : 1 
+ NU*( (U) :1) :1 
+ NU/R/R*( (U) :2) ;2 

- NU/R/R*(V):2 

- NU/R/R/R*(R*V) :2 

- NU/R/R*U 
NU/R*(U):1 


LENGTH VELOCITY 

U* (R*V) :1 0 2 

+ 1/R*V*(R*V) :2 0 2 

= - ( P ) : 2 0 2 

+ NU*( (R*V) : 1 ) ;1 0 2 

+ NLi/R/R*( (R*V) : 2) : 2 0 2 

- NU*TV) ;1 0 2 

+ NU/R/R*(R*U> :2 0 2 

+ NU/R* (U) :2 0 2 

- NU/R*V n 2 



Figure 4. The Navier-Stokes equations In polar coordinates 



In addition to the equation, TATTR was supplied with attribute 
information similar to that for the run of Figure 4 plus data 
on an order of magnitude attribute, as follows: 

Length Velocity Delta 

U 0 1 

W 0 1 

P 0 2 

NU 11 

R 10 

:1 -1 0 

:3 -1 0 

To understand the last column it is easiest to reinterpret the 
equation as nondimenslonal, NU standing for the reciprocal of 
the Reynolds number. (It seems strange at first to assign 
dimensions to the factors and at the same time to consider them 
nondimenslonal, but it is perfectly legal for both interpre- 
tations to coexist.) It is recognized that in a jet at high 
Reynolds number the radial component of velocity is small com- 
pared to the axial. Further, the radial extent of the flow due 
to the jet is small compared to its axial extent, so that changes 
(derivatives) in the radial direction are large compared to 
changes in the axial direction. Elementary boundary layer theory 
teaches us that if the reciprocal of the Reynolds number is of 
order 6^ where 6 is small, then the radial velocity and 
extent are of order S compared to the corresponding axial 
quantities. The last column above, then, gives the exponents of 
6 that indicate the order of the various quantities. 

The result of the TATTR run is shown in Figure 5; is 
seen that the terms are dimensionally consistent but that one is 
small (of order 6^ ) compared to the others. In a succeeding 


1 

0 

0 

2 

1 

-1 

0 
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LENGTH VELOCITY DELTA 

U*(W):1 
+ W* ( W ) ; 3 
= - ( P ) : 3 
+ NU*( ( W) : 1 ) : 1 
+ NU*((W);3):3 
+ NU/R* ( W) : 1 


1 TERMS DROPPED DUE TO DELTA CRITERIA/. LEAVING 5 


LENGTH VELOCITY DELTA 

U* (W) : 1 

+ W* ( W ) : 3 

= -(P) :3 -1 

+ NU*( ( W) : 1 ) : 1 -1 

+ NU/R* ( W) : T -1 


Figure 5. The axial momentum equation In cylindrical coordinates 
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2 0 
2 0 
2 0 
2 0 
2 0 



(optional) stage TATTR dropped the small term. If the relative 
orders of the factors were not known, the run could have been 
set up with separate ’’deltas" for various factors (TATTR can 
handle up to sixteen attributes at once). Various linear 
combinations could then have been tried as criteria for dropping 
terms. With a modicum of physics guiding the combinations to 
be used, this approach would have led to the same results. 

Again, only a hint of the usefulness of TATTR is possible 
here. For the problem considered in Reference 6 the 9843 terms 
(in l8 equations) produced by TENSR were reduced to 529 terms 
by TATTR. These statistics are somewhat misleading — if the 
job had been attempted by hand, most of the negligible terms 
would have been easily spotted. Only the last few hundred 
would have been difficult, but that's plenty to Justify the use 
of TATTR. 


5 . OUTLOOK 


A.R.A.P, is in the midst of switching computers. Our old 
one, a DSC META-4, uses, basically, IBM II 30 software. Our 
new one, a PDF 11/70, uses DEC software. The IBM software is 
geared to the EBCDIC code for character handling, the DEC soft- 
ware to the ASCII code. For TENSR and TATTR, which mainly 
manipulate strings of characters, the difference is significant; 
the conversion will be somewhat more laborious than for ordinary 
A.R.A.P. programs such as the ones that have been developed 
with the help of TENSR and TATTR. 

However, while converting, we will use the opportunity to 

make several Improvements. The principal ones contemplated 
are : 

(1) TENSR and TATTR will be combined. This will allow 
greater convenience for the user and greater efficiency 
of operation. 

(2) The specification of different ranges on different 
dummies will be made possible. 

(3) Conversion to physical coordinates, if desired, will 

be made automatic. In this case, specification of 
symbol substitutions will only be needed in terms of 
physical components Instead of in terms of contra- 
variant components, covariant components, and combi- 
nations thereof . ^ ^ 

(4) The operation of the system will be speeded up and 
thus made more convenient. This will occur partly 
because of the switch to a faster machine, but also 
because of the Implementation of ways of making it 
more efficient which have been recognized since the 
last major revision. 

■ '- 2.9 ■■■■" 


An additional improvement would be to allow substitutions 
before expansion. For example, given 

C'lJK = A*'L»IJ*B'LK + D»IK,J 

and another expression where, say, C”JL’L,K occurs, the 
right-hand side of the former could be substituted in the 
latter with proper shuffling of the indices. This addition 
is quite feasible but is more involved than the other changes 
mentioned. For one thing, a means of handling conflicts in 
indices, as in the example Just given, must be devised. 

Another possibility is to provide output using true super- 
scripts and subscripts. This would depend on the availability 
of suitable output devices and proper software. 

Although TENSR and TATTR were developed largely with NASA 
and Air Force support, they were intended to help in developing 
programs. In some cases these programs were delivered to the 
Government; in other cases the programs were run by A. R. A. P. 
and the results delivered. In any case, TENSR and TATTR were 
not designed for export. The fact that they can be moved to a 
different machine is about to be proved, but whether portability 
will be improved in the process remains to be seen. 

The usage of the combined system is not expected to become 
high in comparison with more general-purpose software, or in 
comparison with programs that churn out numbers. Therefore, 
it is doubtful that the effort to make it portable and, in 
particular, to prepare documentation sufficient to guide a 
remote user would be worthwhile. Instead, I would suggest that 
potential users communicate with A.R.A. P. concerning their 
needs. 



As far as I know, no one else has attempted anything like 
TENSR and TATTR. Coordinate transformations to aid in flow- 
field computations, with varying degrees of automation, are 
plentiful. Reference 1, already mentioned, and Reference 7, 
which proposes conformal transformations in the cross-flow 
plane for computing supersonic flows, are just two recent 
examples of many papers describing many approaches. 

Pormac and similar systems do algebra and calculus auto- 
matically— and are much better at it than TENSR and TATTR — 
but have no provision for incorporating tensor notation. 

Howard (References 8 and 9» among other papers) has used such 
systems for doing tensor analysis, but he must write out 
covariant derivatives in terms of Christoffel symbols (for 
example, the right side of equation (8) instead of the left), 
and he uses explicit nested DO loops for each pair of dummy 
indices. But he recognizes the need for making such processes 
automatic, not only for conserving manpower, but also for 
drastically reducing the likelihood of making mistakes. 
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6. NOTE FOR TENSOR EXPERTS 


This section covers some points that may be of Interest 
to those with a relatively thorough knowledge of tensors. 

First, TENSR Is misnamed. It operates on the symbols 
supplied to It, as described In Section 3» whether they 
represent tensors or not. In fact, as soon as It Introduces 
the Chrlstoffel symbols It operates on the various terms 
separately, and the separate terms are not tensors. In other 
words, TENSR Is not restricted to tensors. 

Second, there Is no need for the space Involved to have 
a metric. The metric Is relevant only If It appears In the 
Input equations, either explicitly or Implicitly through the 
use of ! . This Is one reason the Chrlstoffel symbols (more 
properly Interpreted as affine connections when there Is no 
metric) are separate Inputs to TENSR. (The other reason Is 
that to be able to determine them from a metric, TENSR would 
have to be taught to differentiate.) 

Third, expressions can be expanded without anything being 
known about the space except Its dimensions and that It has 
affine connections and. If used, a metric (but not the values 
of the affine connections or the metric), as is illustrated 
in Figure 2. Further, it Is possible to suppress the intro- 
duction of the affine connections (Chrlstoffel symbols), or 
to suppress the expansion of dummy indices, or both. Thus, 
there is great flexibility in determining what is to be done. 

Fourth, the number of dimensions is limited to ten solely 
because that Is the number of numerical digits available In 
computer output devices. The user may select any set out of 
the ten for dummy expahsions and for free indices. In particular 
for runs Involving four-dlmenslonal space-time, the user is free 
to let whim guide the choice rather than being forced to use 
0,1, 2,3 or 1,2, 3,4 or even being forced to choose between 
them. ■ ■ 


7. bibliography 


This section contains a few suggestions for those who 
wish to pursue the subject of tensors, I have made no effort 
to Inspect everything that is available; I report only on some 
works that I have chanced upon or that have been brought to 
my attention by others. 

The first group considers tensors as tools useful in 
various branches of applied mathematics; if they consider 
general relativity, it Is as one of several applications. 

The rest are all devoted to general relativity and develop 
the theory of tensors for that application only. 

Mlchal, Reference 10, was a pioneer in making tensors 
accessible to engineers, but Sokolnikoff, Reference 4, is more 
thorough both in the basics and in the applications. A modern 
treatment by Arls, Reference 5 , is very readable; the appli- 
cations are restricted, as the title indicates, to fluid 
mechanics. Budiansky provides a remarkably thorough treatment 
in less than 50 pages (Reference IT). He includes Just a 
taste of applications by considering mechanical deformations. 

It is no wonder that much of what is written on tensors 
is in works on gravitation. Although they are extremely useful 
in many other fields, tensors are deemed indispensable in 
general relativity. Consider first three short works by men 
whose names have become bywords. Foremost among these is the 
monograph. Reference 12, by Einstein. Pear not! The parts 
dealing with mathematics are clear and straightforward, (I’ll 
make no comment on the physics. ) Schrodinger, Reference 13 , 
is dry and elegant with barely a; hint of what it ’ s all about . 

By contrast, Dirac is discussing black holes by page 32 of 
Reference l4. 
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Two more thorough treatises worth looking into are 
Anderson, Reference 15» e.nd Weinberg, Reference 16, 

Finally, there is the magnificent Ovavitation (Reference 
17 ) by Mlsner, Thorne, and Wheeler, This book has exhaustive 
treatments of a multitude of subjects (the discussion of 
black holes doesn’t start until page 872), it Includes, of 
course, a thorough study of tensors, but also includes a 
parallel development of the mathematics from the viewpoint 
of differential geometry, along with discussions of the 
relationships between these two subjects. In spite of its 
bulk and comprehensiveness it is lively and accessible. 
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